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Collisional losses, decoherence, and frequency shifts 
in optical lattice clocks with bosons 
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Phystkalisch-Technische Bundesanstalt, Bundesallee 100, 38116 Braunschweig, Germany 
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We have quantified collisional losses, decoherence and the collision shift in a one-dimensional 
optical lattice clock with bosonic 88 Sr. The lattice clock is referenced to the highly forbidden 
transition 1 So— 3 Po at 698 nm, which becomes weakly allowed due to state mixing in a homogeneous 
magnetic field. We were able to quantify three decoherence coefficients, which are due to dephasing 
collisions, inelastic collisions between atoms in the upper and lower clock state, and atoms in the 
upper clock state only. Based on the measured coefficients, we determine the operation parameters 
at which a fD-lattice clock with 88 Sr shows no degradation due to collisions on the relative accuracy 
level of 10 -16 . 

PACS numbers: 06.30.Ft, 32.70.Jz, 34.50.Cx, 37.10.Jk 
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Laser spectroscopy of transitions with millihertz 
linewidth has become the essential tool for the rapidly 
evolving optical clocks based on single ions [l[ or neutral 
atom clouds These clocks promise a realization of 
the base unit 'second' with higher stability and accuracy 
than the currently employed frequency standards using 
a microwave transition in 133 Cs. Optical clocks further- 
more enable tests of new theories, which predict e.g. the 
variation of fundamental constants Q, or they can be 
used for relativistic geodesy field experiments Q. Due 
to long coherence times and very sensitive interrogation, 
optical clocks with neutral atoms are also ideally suited 
to explore the physics of ultra-cold collisions. The inter- 
pretation of collision experiments provides insight into 
the coupling between the atoms in different electronic 
states at ultra-low temperature. Understanding collision 
properties of alkaline earth atoms is required for proposed 
quantum computing schemes [B| and to test atomic struc- 
ture calculations p. 

Strontium is widely used for magic wavelength optical 
lattice clocks with, both, bosons and fermions [3, [D, [8|, [9| • 
At ultra low temperature, collisions of fermions are 
strongly suppressed by the Pauli exclusion principle [Io| . 
The bosonic isotope 88 Sr offers higher natural abundance 
and easier laser cooling. It is therefore of interest espe- 
cially for transportable setups. Collisions are suppressed 
if a three dimensional lattice with an occupancy of at 
most one atom per site [8| is used but the complexity of 
setup and experimental sequence is increased. The use of 
88 Sr in a ID-lattice clock can therefore be advantageous 
if the effects of collisions were well quantified 0, Hlf . 

Three effects induced by collisions are expected: in- 
elastic losses, collision broadening of the transition, and 
a collision-induced frequency shift. We have quantified 
all three in our lattice clock setup. With our results we 
are able to predict operation conditions of a 88 Sr lattice 
clock under which collisions are not expected to limit the 
relative accuracy on the level of 10~ 16 . Additional shifts 
due to clock-laser induced ac Stark effect and quadratic 



Zeeman effect in bosonic lattice clocks are theoretically 
well understood and can already be controlled with high 
precision. For these, a relative uncertainty far better 
than 10 -16 is expected for the near future [121 ]. 

Our experiment was described before [l3j ]. In short, 
we first prepare the 88 Sr atoms in a magneto-optical trap 
(MOT) operated on the strong 461 nm transition 1 Sq — 
1 Pi. The atoms are loaded for up to 250 ms from a 
Zeeman-slowed atomic beam, which is deflected into the 
MOT region by means of an optical molasses. 

The atoms are further cooled in a second cooling phase 
on the 689 nm intercombination line 1 So — 3 Pi- First, 
the cooling light is frequency modulated to enlarge the 
capture velocity range of the second stage MOT. Then, 
after 70 ms, the modulation is turned off and within an- 
other 70 ms the atoms are further cooled to below 4 /jK. 
A horizontal ID-optical lattice is overlapped during the 
whole cooling cycle with the MOT. It is created by a laser 
beam of up to 300 mW power at the magic wavelength 
of 813 nm. The beam is focused to a waist of 34 /jm and 
retro-reflected into itself. The waist size is confirmed via 
measurements of trap oscillation frequencies. After the 
second cooling phase, up to 3 x 10 6 atoms are trapped 
in the optical lattice (about 1000 atoms per site). The 
atom number and temperature were determined by ab- 
sorption image measurements after a variable expansion 
time. Untrapped atoms are falling out of the interroga- 
tion area within 150 ms. Then the trapped atoms are 
interrogated on the clock transition 1 So— 3 Po at 698 nm. 
The clock laser beam of 2.5 mW maximum power is over- 
lapped with the lattice through a dichroic mirror and fo- 
cused to a waist radius of about 39 /im at the position 
of the atoms. To induce a dipole transition matrix ele- 
ment on the clock transition, we apply a homogeneous 
magnetic field of up to 3 mT. 

After excitation on the clock transition the remain- 
ing ground state atoms 1 So are detected by recording 
the fluorescence during a 461 nm MOT phase and then 
pushed away by resonant 461 nm radiation. In the same 
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FIG. 1: (color online) Inelastic losses by collisions between 
3 Po atoms (open circles) and a mixture of 3 Po and 1 So atoms 
(full symbols). Points are the measured data, lines are the 
results of fits. The inset shows schematically the involved 
atomic levels. The initial number of ground state atoms is 
about 3 x 10 6 . 



way, excited 3 Po atoms are detected after optically pump- 
ing them to the 3 Si state from where they decay to the 
ground state via the 3 Pi state (inset Fig. [TJ). Photoas- 
sociative losses by the repumping light are avoided by 
switching off the lattice and expanding the cloud of 3 Po 
atoms before repumping 

We first address inelastic collisions: Without excita- 
tion on the clock transition, we observe a background 
gas limited lifetime of the atoms in the lattice of 7 s to 
8 s. In the presence of atoms in the 3 Po state, density 
dependent shortening of the trap lifetime is observed. In- 
clastic two-body losses are possible by collisions of 3 Po 
atoms with either 3 Po or 1 So atoms. The decay can be 
modeled by the coupled differential equations 



Pc = -rp c - IgePgPe 
Pg = -I>g - 7gcPePg 



7ceP c 



(1) 



Here, p g (p c ) denotes the local atomic density in the 
ground (excited) state, T is the inverse trap lifetime and 
7gc (7ce) the loss coefficient for ground - excited (excited 
- excited) state collisions. Radiative decay is much slower 
and neglected. 

To distinguish experimentally between both inelastic 
collision channels, we first removed the ground state 
atoms from the lattice and observed the decay of the 
3 Po population only (Fig. [1]). In this case, the Eqs. [TJ 
simplify, decouple, and can be solved analytically. Inte- 
gration over the spatial coordinates of a single lattice site 
provides an expression for its atom number: 



N(t) = N - 



exp (-Tt) 



1 + iVo7ee/(7r 3 / 2 r W > z ) [1 - exp(-rt)] ' (2) 

Nq denotes the initial atom number in a lattice site, w r 
and w z are the 1 / e 2 -radii of the atomic cloud in a lattice 
site at a given temperature and set of lattice parameters. 



In our experiment many independent lattice sites are 
occupied with different atom numbers. We numeri- 
cally sum Eq. [5] over the lattice sites. The initial ax- 
ial population distribution is described by a Gaussian 
of 1/e 2 radius w' (determined from absorption images). 
From several data sets we determined the loss coefficient 
7c = (2.0 ± 0.7) x 10~ 18 m 3 /s. This is in good agree- 
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ment with the value given by Traverso et al. 
uncertainty is dominated by the uncertainty of the calcu- 
lation of the radii Wi and the scatter of the fit parameters. 
We neglect evaporation losses, since the lattice depth of 
&b ■ 45 pK is much larger than the atomic thermal energy 
of below /cb • 4 pK. 

We then investigated a mixed sample of 3 Po and 1 So 
atoms. Collisions 3 Po+ 1 So are responsible for the dif- 
ferent decays visible in Fig. [TJ We determined the loss 
coefficient 7 go = (2.3± 1.1) x 10~ 19 m 3 /s by numerical in- 
tegration of Eqs. [TJ Contributions to the uncertainty are 
again the determination of the W{ and details of the fit. 
Careful consideration of the fitting strategy is required 
because losses are dominated by 3 P + 3 Po collisions as 
analyzed above. Best accordance of data and fit was ob- 
tained by fitting 7 ee together with 7 ge . Nevertheless, we 
constrained 7 ee to the value obtained under the clean con- 
ditions of the first set of measurements. We also found 
slightly different results for 7 ge depending on whether we 
fitted only the excited state population, the ground state 
atom number or both simultaneously. We attribute this 
scatter to inelastic collisions that diminish the 3 Po atom 
number during the detection of the ground state atoms. 
We corrected the excited state population for collisions 
with 3 Po atoms but neglect other types of collisions. The 
scatter of 7 ge determined by different strategies is finally 
represented by the uncertainty given above. 

Different loss mechanisms are possible at the respective 
collision asymptotes. The asymptote 3 Po+ 3 Po correlates 
at short internuclear distances to a state, which can 
decay by fluorescence to lower molecular states. More 
likely as source for the losses are non-adiabatic transi- 
tions to molecular potentials dissociating to lower asymp- 
totes. In the case of collisions at 1 So+ 3 Po, the latter 
path seems to be very unlikely. Only one potential disso- 
ciates to the one lower asymptote 1 So+ 1 So and it is very 
steep at the energy of the 1 So+ 3 Po asymptote. A direct 
spontaneous decay to that ground state potential is not 
possible due to selection rules but at short internuclear 
distance complex coupling exists to molecular states that 
are subject to spontaneous decay 
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rect pathway can be seen as a qualitative explanation for 
the significantly smaller loss coefficient at the asymptote 
1 So+ 3 Po compared to the asymptote 3 Po+ 3 Po- How- 
ever, we cannot exclude three-body losses or processes 
involving a lattice photon. 

As second decoherence process in a 88 Sr lattice clock 
we identified strong density dependent broadening of the 
clock transition (see Fig. [2J) . To model the observed spec- 



3 




clock laser detuning (kHz) 



clock laser detuning (kHz) 



FIG. 2: (color online) Rabi oscillations (top row) and spectra 
of the clock transition (below) for different atom numbers 
indicated on the ordinate. The solid lines are calculated with 
the density matrix model. The red curves stem from fits of the 
dephasing constant 7d cp , while the blue ones are a calculations 
with 7dc P = 0. 



tra we extended the differential equations Eqs.Q]to a for- 
malism of the full density matrix p (pu — p g , P21 = p c ) 
including the excitation dynamics, and integrated them 
numerically: 



p = --[H,p]+n(p). 



(3) 



In the rotating wave approximation, the atom-light 
Hamiltonian and the relaxation matrix IZ(p) become 

H / o n/2 \ 
h ~ \n/2 a ) 

R(p)n = -(r + 7gcP22)pn 

Tl(p)l2 = K(p)n = ~[7ge(Pll +P22)/2+7cc/022/2 

+r + L + 7do P /Oii]pi2 

fc(p)z2 = ~(r + 7gePll + 7ccP22)P22- (4) 

The Rabi frequency f2 can be calculated from the ex- 
perimental parameters 0, [l~5| , the detuning A is deter- 
mined from the spectrum itself. Dephasing or damping 
of the coherences in the density matrix is introduced by 
the coefficient L representing an effective laser linewidth 
and inhomogeneous damping due to spatial dependent 
Rabi frequencies. The dephasing by elastic collisions is 
described by 7dep- We assume dephasing by collisions 



with ground state atoms. Different mechanisms can be 
imagined but the experimental data do not allow to dis- 
tinguish between them. 

For the sake of fast computation, we neglected spatial 
variations in the excitation probability due to the lattice 
or the clock laser profile. Furthermore, a rectangular 
population distribution of full width 2w' in axial direction 
was assumed. We have verified in the case of the inelastic 
losses that the effect of this approximation is minor. 

To fix the excitation dynamics we fitted observed Rabi 
oscillations and a scan of the clock transition at an atom 
density, which does not significantly broaden the tran- 
sition. The fitted Rabi frequency is about 17% smaller 
than the one calculated from the peak clock laser in- 
tensity and the homogeneous magnetic field, an effect 
we attribute to the finite Lamb-Dicke parameter. The 
FWHM laser linewidth was fitted to be about 36 Hz, 
which is roughly consistent with a transition linewidth of 
35 Hz we have observed at an excitation pulse length of 
35 ms and a few thousand atoms. To match the obser- 
vations at higher density we had to include (besides the 
inelastic collisions determined before) a dephasing coef- 
ficient of 7dc P = (3.2 ± 1.0) x 10~ 16 m 3 /s (see Fig. EJ. 
The effect is stronger by a larger factor (10 3 ) compared 
to collisions 1 S + 3 Pi in a 1000 K Sr gas [lj|, which could 
indicate the existence of a scattering resonance. We did 
not observe a significant temperature dependence of any 
collision coefficient 7 between 2 pK and 4 pK. 

The third collision effect being very important for the 
accuracy of an optical clock with neutral atoms is the 
frequency shift 5v due to collisions. We have measured 
5v by two stabilizations with interleaved probing cycles 
[13, EH- The clock laser probes either about 10 4 atoms 
(used as reference) or a higher atom number (and thus 
density). This is systematically varied via the duration 
of the initial MOT stage. The density shift is revealed 
as difference of the offset frequencies for the alternated 
stabilizations between the clock laser and its reference 
cavity. We have verified that no systematic frequency 
shifts occur for equal atom numbers on the level of a few 
10 mHz. We shortened the cycle time to r cyc i c = 200 ms 
and detected the excited state population only. This way 
the stability of the difference relative to the optical fre- 
quency improved to be better than 10 -15 at an averaging 



time r = 10 s and scaled as l/yr. The dependence of 
the observed frequency shift as a function of the differ- 
ence of atom number in the two stabilizations is shown 
in Fig. [3l At low density we observe a linear depen- 
dence of the shift on the atom density; for high density 
we observe saturation of the shift. We attribute this to 
the changed dynamics during the excitation due to losses 
and dephasing. A theoretical analysis of the excitation 
process including a excitation dependent frequency shift 
will be carried out in the future. The clock laser was 
stabilized by probing the atoms with a 5 ms 7r-pulse de- 
tuned such that the excitation probability was reduced 
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FIG. 3: (color online) Observed density dependent shift of 
the 88 Sr clock transition. The open symbols were excluded 
from the linear fit (full line, see text). The inset shows the 
dependence of the shift on the excitation probability with a 
linear fit. The triangle in the main graph represents the data 
in the inset. The shift was rescaled according to the trap 
parameters and atom temperature. 

to about 35% after the pulse. 

From a linear fit through the origin (excluding the open 
data points in Fig. [3] because of saturation), the coeffi- 
cient of the absolute collision shift as a function of the 
atom number iV can be determined to better than 4%. 
The coefficient for the shift as function of the density is 
Av p = 2Swk 2 w'w^w z /NX = (7.2 ± 2.0) x 10" 17 Hz-m 3 
with an uncertainty mainly due to the density determina- 
tion (A: magic wavelength). A temperature dependence 
of the shift was not observed between 2 /iK and 4 /iK. 

The inset in Fig.[3]shows the dependence of the density 
shift on the excitation probability after the clock laser 
pulse, which was varied by the clock laser detuning in the 
stabilization cycle. We observe a weak relative influence 
of (2.9 ± 4.5) x 10~ 3 per percent excitation probability. 

Having quantified three collision influences, we give 
guidelines for the design of a ID-lattice clock with 
bosonic 88 Sr. Assuming a typical lattice depth of 
/cbTO /iK, an atom temperature of 3 /iK, and an avail- 
able lattice laser power of 300 mW, one could choose a 
lattice waist of 75 /tm. With w' = 280 /im and at the 
current level of accuracy for the density shift correction 
of 4%, a density shift of about 1 Hz would be tolerable to 
reach a relative accuracy of 10~ 16 . This would limit the 
atom number to about 2 x 10 4 , a value comparable to or 
larger then in present lattice clocks with 87 Sr. The col- 
lisional broadening is then about 1.3 Hz. A new density 
shift measurement in the proposed lattice should yield an 
improved correction and allows for increasing the atom 
number until the collisional broadening becomes relevant. 
Aiming at a line width of about 10 Hz, operation with 
more than 10 5 atoms is feasible. With a cycle time of 
200 ms, the quantum projection noise limited relative 



stability in one second is 2 x 10~ 17 . At this density, ac- 
cording to our findings losses do not distort the observed 
line or limit the excitation probability. 

We conclude that collisions do not impose limitations 
for a bosonic 88 Sr lattice clock compared to state-of-the 
art ones with 87 Sr. We expect that these newly observed 
and quantified decoherence mechanisms also will have im- 
pact on the design and operation of experiments with al- 
kaline earth atoms that require long coherence times like 
quantum computing, millihertz linewidth lasers (l9j . or 
optical Feshbach resonances |20| . 
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